The Large iV Reduction in Matrix Quantum Mechanics 
— a Bridge between BFSS and IKKT — 
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The large N reduction is an equivalence between large N gauge theories and matrix models dis- 
covered by Eguchi and Kawai in the early 80s. In particular the continuum version of the quenched 
Eguchi-Kawai model may be useful in studying supersymmetric and/or chiral gauge theories non- 
perturbatively. We apply this idea to matrix quantum mechanics, which is relevant, for instance, 
to nonperturbative studies of the BFSS Matrix Theory, a conjectured nonperturbative definition of 
M-theory. In the bosonic case we present Monte Carlo results confirming the equivalence directly, 
and discuss a possible explanation based on the Schwinger-Dyson equations. In the supersymmetric 
case we argue that the equivalence holds as well although some care should be taken if the rotational 
symmetry is spontaneously broken. This equivalence provides an explicit relation between the BFSS 
model and the IKKT model, which may be used to translate results in one model to the other. 
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PACS numbers: ll.25.-w; 11.25.Sq 

Introduction. — The large N reduction states that 
large N gauge theories are equivalent to matrix models 
obtained by dimensionally reducing those theories to a 
point. Recent developments in the exact calculation of 
the low-energy superpotential in J\f = 1 supersymmetric 
gauge theories |2| may also be regarded as a particu- 
lar application of the large N reduction Q. We apply 
the original idea to matrix quantum mechanics, which is 
relevant, for instance, to nonperturbative studies of the 
BFSS Matrix Theory Q proposed as a nonperturbative 
definition of M-theory. The same model may also be in- 
terpreted as a low-energy effective theory of DO-branes in 
type IIA superstring theory. (In the cases where fermions 
are absent, interesting results are already obtained by a 
direct lattice approach [jj. See Ref. for a lattice for- 
mulation in the supersymmetric case.) 

In Ref. 0] it was conjectured that the U(iV) gauge 
theory on an infinite lattice is equivalent in the large N 
limit to a one-site model called the Eguchi-Kawai (EK) 
model 0. This conjecture was based on the observation 
that the Schwinger-Dyson equations of the two theories 
coincide under the assumption that the U(l) d symme- 
try Ufj, — - > e lafi V 'fj, of the EK model remains unbroken. 
The U(l) d symmetry, however, turned out to be sponta- 
neously broken in the weak coupling region 0, , and the 
quenched EK (QEK) model was proposed as a remedy 
[3. The "quenched" eigenvalues of were interpreted 
later as lattice momenta 9j, and this interpretation en- 
abled an extension to non-gauge theories. The equiva- 
lence for the QEK model was confirmed both perturba- 
tively an d nonperturbatively 0, El ■ 

The continuum version of the QEK (cQEK) model has 
been proposed in Ref. 0] . In spite of its potential use- 
fulness in studying supersymmetric and/or chiral gauge 
theories, the EK equivalence for the cQEK model has not 



been checked directly so far. 

The twisted EK model 0] , which was proposed as an 
alternative way to cure the problem of the original EK 
model, has been interpreted recently as field theories on 
a noncommutative geometry |l4| . However, this model 
is not useful for our purpose since the twisting proce- 
dure is applicable only to even space-time dimensions, 
and moreover the continuum version can be defined only 
in the strict N — 00 limit 15]. As another recent de- 
velopment, it has been found that the breaking of the 
U(l) d symmetry can be avoided by keeping the size of 
the lattice sufficiently large, and that the critical size in 
physical units becomes finite in the continuum limit |l6| . 

In this letter we consider the cQEK model for ma- 
trix quantum mechanics, and discuss its EK equivalence. 
In fact the corresponding cQEK model is analogous to 
the IKKT model |l7j proposed as a nonperturbative for- 
mulation of type IIB superstring theory. Based on this 
observation, we discuss an explicit relation between the 
BFSS model and the IKKT model. 

The models. — We study a one-dimensional U(iV) 
gauge theory, or matrix quantum mechanics (QM), with 
adjoint scalars Xi(t) (i = 1, • • • , d — 1), which are N x N 
Hermitian matrices. The action is given by 
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where DXi — J^-X^ — i [A, Xi\ represents the covariant 
derivative with A being the Id gauge field. The model 
can be obtained formally by the dimensional reduction 
of (i-dimensional Yang-Mills theory to one dimension. In 
what follows we fix the gauge to A(t) = 0. This model is 
UV finite and hence it does not require any UV regular- 
ization. As a result, the parameter g in Q can always 
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be scaled out by an appropriate rescaling of the matrices 
and the time coordinate t. We take g = which turns 
out to be convenient when we discuss the large N limit. 

Applying the quenching prescription 0, Il2| to the ma- 
trix QM QJ, we obtain the corresponding cQEK model 

<S C QEK,f2 = ~Netr Q [r»,X ( ] 2 + j [JC<,X,] 2 ) , (2) 

where Xi (i = 1, ■ ■ ■ , d— 1) are NxN traceless Hermitian 
matrices and Q, = diagfwi, • • • ,lon) represents frequency 
variables to be integrated over the region [— 4, 4j, where 
A = — . (More precisely, the difference oji — ujj corre- 
sponds to the frequency.) 

As U(A) invariant operators in the matrix QM QJ, 
let us consider O = j^tv X il (ti) ■ ■ ■ X in (t n ) and the cor- 
responding operator O = -^tr X^ti) ■ ■ ■ X in (t n ) in the 
cQEK model Q, where X t {t) = e int X ie ~ int . Then, 
according to the EK equivalence, 

lim (0) m QM = lim lim / dtt (0) C QEK.n , (3) 

N—>oc A^oo N— >oo J 

where ( • } m QM and ( • ) C QEK,n represent the expectation 
values with respect to the actions and J5J), respec- 
tively, and / dfl = njli Ha The order of the two 
limits on the r.h.s. of © cannot be inverted. 

Monte Carlo simulation. — We test the equivalence 
(0 directly by Monte Carlo simulation. In order to sim- 
ulate the Id model (Q, we discretize the time "t" , and 
obtain the action 

\ n—l ^ ' 

-^ritrlX.H,^)] 2 ] , (4) 

n— 1 ^ 

where a represents the lattice spacing. The NxN Her- 
mitian matrices Xi(n) (n — 1, • • ■ , T) sit on each site, 
and we do not impose periodic boundary conditions. In 
order to retrieve the original model we have to take 
the a — > and r = a T — > oo limits, where the order of 
the limits should not matter. The system is invariant un- 
der Xi(n) <— > Xi(n) + Pi 1, and we fix the corresponding 
zero mode by imposing the condition * trXj(n) = 0. 

Monte Carlo simulation of the models 10} and @ can 
be performed by the heat bath algorithm [18|. For the 
cQEK model, we first take an average over the ensemble 
{Xi} with the action J5J for each Q, and then perform the 
averaging over f2 by using random numbers distributed 
uniformly within the region [— 4>t]- Correspondingly 
there are two types of statistical errors. We estimate the 
total error by taking an average over the error of the first 
type and adding the one of the second type. 



We focus on the d = 4 case in what follows. Let us 
first consider the observable 0\ — j^ELi' r ft(' 1 ) 2 ) 
for the lattice model |@}. In Fig. ^ we plot the results for 
(Oi) lat against a for t = 5.0, 7.0, 10.0, 15.0 at N = 16. 
We fit the results to (Oi)i a t = C\ + C 2 a + C 3 a 2 , from 
which we obtain the continuum limit. 
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FIG. 1: The observable (Oi)i a t for the lattice model (QJ with 
N = 16 is plotted against a for r = 5.0, 7.0, 10.0, 15.0. The 
lines represent fits to (Ci)i a t = C\ + C2 a + C3 a 2 . 
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FIG. 2: The observable (Oi)iat for the lattice model (0J after 
taking the a — -> and r — > 00 limits are plotted against ^3-. 
The straight line represents a fit to Ce + C7 A r_2 . 

The results obtained in the continuum limit for var- 
ious t can be fitted to lim a _o(Oi)iat = C4 + C^t^ 1 , 
from which we obtain the r — > 00 limit. We redo this 
analysis for N — 20, 32. The results can be fitted to 
Cq + Cj N~ 2 as one can see from Fig. [5] This large N 
behavior is analogous to the one observed in the bosonic 
IKKT model [18j • Thus we can obtain the final result for 
the matrix QM 0J. 

Let us move on to the cQEK model @. The corre- 
sponding observable is <D\ = -^-tr (A 2 ). In Fig. [3]we plot 
the results (Oi) c qek obtained for N — 32, 48, 64 with 
e = 0.2, 0.25, 0.3, 0.5. For fixed e the large N behavior is 
given by Cs + Cg N~ 2 . The coefficient Cg, however, di- 
verges as -4 ■ This can be understood if we recall that the 
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density of in the region [— — , |-] is p = ^ = Since 
finite TV effects at small e come mainly from the finiteness 
of p, any power of N should be associated with the same 
power of e. This motivated us to plot the results against 
jj^rp in Fig- 01 where the results indeed lie on a single 
straight line for e < 0.3. Fitting all the data for e < 0.3, 
we obtain the final result for the cQEK model (0). 



observable 


matrix QM 


cQEK 


<^trX^ 

( lv tr ^ "^J ) 
{^trpGX,) 2 ) 


1.279(1) 
2.15(1) 
1.23(2) 


1.28(1) 
2.11(5) 
1.27(4) 



TABLE I: Comparison of the Monte Carlo results obtained 
for the matrix QM and the cQEK model @. 
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FIG. 3: The observable (Oi> c qek for the cQEK model is 
plotted against nvlp ^ or 6 = 0-25, 0.3, 0.5. The straight 

line represents a fit to (Oi) c qek = Cg + Cio (Ne)~ 2 . The 
cross on the vertical axis (the error bar is almost invisible) 
represents the final result for the matrix QM (0. 

We repeat this analysis for other observables and sum- 
marize the results in Table U where we find good agree- 
ment within error bars. We consider this as a compelling 
evidence for the EK equivalence |J5}. 

Schwinger-Dyson equations. — Let us discuss a pos- 
sible explanation for the equivalence O) based on the 
Schwinger-Dyson (SD) equations [lL [T3. Il9j| . 

In the matrix QM we consider 

(tr (X a X il {ti) ■ ■ ■ Xi n (t„))) m QM , 

where A a represents a generator of the U(N) group. 
Changing the integration variables Xi(t) i— > Xi(t) + 
X a Si t i S(t — to), summing over a, using the identity 
S Q (^ t %( Aa )fc' = faSjk, and assuming the large TV 
factorization property (OO')mQM ^ (0)mQM<O')mQM, 
we obtain a closed set of equations for the observables 
(C)mQMi which is the SD equations. 
In the cQEK model we consider 

J dn(ti(e int n a e- im °X n (ti)---X ln (tn)))cQEK,n ■ 

Changing the integration variables X{ i— > Xi + ^X a 6i : i , 
we proceed similarly to the matrix QM case assuming in 
addition that the ^-integration factorizes as well. In this 
way we obtain a set of equations, which are the same as 
the SD equations for the matrix QM with the identifica- 



tion © except for the terms (s = 1, ■ • ■ , n) 
U io>i i rfSl(tr {e tn{t °-^X n (h) ■ ■ ■ AVi (* s -i)))cQEK,n 
x y"dr!(tr(A > l3+1 (i s+1 )---A > l „(t n )e- lO(t =- t0) )) C QEK,o , 

whose counterpart exists in the matrix QM only when 
t s = to. If we neglect the fact that the integration region 
of SI is bounded by A, the system is invariant under the 
shift SI i— > fl+a 1 except for the factors e ± in (*s-*o) ; which 
make the two SI- integrations vanish separately for t s =/= to- 
Since the model we are studying is UV finite, we may 
expect that the integrand for the Sl-integration is deter- 
mined by the u/j's which are clustered in a finite region. 
Then the integration over the position of the cluster gives 
a factor of f(t s - to) = ^e^^"* ), and we obtain 

the coefficient \f(t s — to) 2 — ► S(t s — to), which agrees 
with the corresponding one in the matrix QM in the 
e —> limit. Assuming that the SD equations have a 
unique solution, we obtain the equality ©. 

Relation to the IKKT model. — The cQEK model © 
for the matrix QM has the action of the same form as 
the bosonic IKKT model H3 Smkkt = -f tr [ A m A A 2 
if we identify A, = e 1 / 4 X l (i < d - 1) and A d = e 1 / 4 ^. 
The only difference is that Ad has to be quenched. In the 
bosonic IKKT model it is known that (^tr (^4 At ) 2 )biKKT 
is of 0(1) 0, and the extent of the eigenvalue distribu- 
tion is of 0(1) in all directions. The quenching stretches 
the eigenvalue distribution of Ad to have an extent of 
0(e -3 ' 4 ), and as a result the eigenvalue distribution of 
Ai (i < d— 1) shrinks to 0(e 1 ^ 4 ) according to our results. 
This implies that, in spite of their formal similarity, the 
bosonic BFSS model and the bosonic IKKT model probe 
totally different regions of the configuration space. 

How about the supersymmetric case? In d = 4, since 
the situation is qualitatively the same (no spontaneous 
breaking of the rotational symmetry and (^tr (A^) 2 ) ~ 
0(1) HlJ), all the statements for the bosonic case should 
equally apply. 

In d = 10 (i.e., the IKKT model), there are certain ev- 
idences that the eigenvalue distribution of A^ collapses 
dynamically to a four-dimensional hypersurface . Let 
us assume further that the eigenvalue distribution ex- 
tends to infinity in four directions at large N, and that 
the distribution is uniform in these directions, as ex- 
pected if the IKKT model really describes our 4d space- 
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time. When we quench the eigenvalues of Aiq, the 4d 
hypersurface will arrange itself to have the required ex- 
tent in the 10th direction. In this situation a vector in 
the 10th direction may have non-zero components in di- 
rections orthogonal to the 4d hypersurface, which implies 
that the system is not invariant under a shift in the 10th 
direction. We therefore suspect that the EK equivalence 
between the BFSS model and the corresponding cQEK 
model does not hold as it stands. 

We may, however, add a "mass term" \N m 2 J dttrX? 
to the BFSS model, and consider the corresponding 
cQEK model, in which the mass term forces the 10th 
direction to be included in the 4d hypersurface. This 
guarantees the translational invariance in the 10th direc- 
tion except at the boundary, which should be fine again 
due to the UV finiteness. Therefore the EK equivalence 
in this case is expected to hold at any finite m, and it 
should hold, too, even if we send m to zero eventually. 

We speculate that actually quenching is not needed 
since the 10th direction, in which we do not have the mass 
term, will extend to infinity as N — > oo anyway. This im- 
plies that the IKKT model with a small mass term in 9 di- 
rections, which should be removed after taking the large 
TV limit, is equivalent (in the sense of Eguchi-Kawai) to 
the BFSS model. In particular if the IKKT model has 4 
extended directions, the BFSS model should have 3 out of 
9 transverse directions extended, thus breaking the SO (9) 
symmetry down to SO (3) spontaneously. This equiva- 
lence is somewhat reminiscent of the T-duality, which 
relates the BFSS model on a circle to the IKKT model 
with Taylor's compactifying condition j^. Note, how- 
ever, that our conjecture is different in that the two mod- 
els can be both considered in the decompactified limit. 

Concluding remarks. — In this letter we have pro- 
vided the first direct confirmation of the EK equivalence 
for the cQEK model in the case of matrix QM or the Id 
gauge theory. Whether it holds also in higher dimensions 
is an interesting open question. This is relevant to non- 
perturbative studies of Matrix String Theory |24|, and 
it is also important in view of the difficulties in formu- 
lating supersymmetric and /or chiral gauge theories on 
the lattice. (See Refs. 0, |3| for recent developments.) 
As is clear from the discussion on the SD equations, the 
question is highlynontrivial in the presence of UV diver- 
gences. In Ref. 12] it is speculated that the equivalence 
holds if the UV divergence is at most logarithmic. A nat- 
ural starting point in this direction would be to study the 
cQEK model for the 2d pure Yang-Mills theory, and see 
whether the known exact results can be reproduced. 
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